Exact wave-packet decoherence dynamics in a discrete spectrum environment 
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We find an exact analytical solution of the reduced density matrix from the Feynman- Vernon 
influence functional theory for a wave packet influenced by an environment containing a few dis- 
crete modes. We obtain two intrinsic energy scales relating to the time scales of the system and the 
environment. Different relationship between these two scales alters the overall form of the solution 
of the system. We also introduce a decoherence measure for a single wave packet which is defined as 
the ratio of Schrodinger uncertainty over the delocalization extension of the wave packet and char- 
acterizes the time-evolution behavior of the off-diagonal reduced density matrix element. We utilize 
the exact solution and the docherence measure to study the wave packet decoherence dynamics. 
We further demonstrate how the dynamical diffusion of the wave packet leads to non-Markovian 
decoherence in such a microscopic environment. 

PACS numbers: 03.65.Yz, 03.65.Db, 42.50.Lc 



I. INTRODUCTION 

Wave packet dynamics has been extensively studied 
since the birth of quantum mechanics Studying de- 
coherence dynamics of a wave packet still attracts atten- 
tion in various research topics, from quantum dissipation 
dynamics, quantum cosmology, quantum measurement 
theory, quantum information processing, to the founda- 
tion of quantum mechanics The process of quan- 
tum decoherence is triggered by the interaction of the 
system with its environment. A widely used model of 
the environment is a set of harmonic oscillators [1, 1^. 
Environments with continuous distribution of oscilla- 
tor frequencies have been extensively studied concern- 
ingthe dissipation effects and the induced decoherence 
M 0, i, i [3, [H [11 M MM,MM,M- We shall 
consider in this paper an environment containing finite 
discrete modes. 

Wave packet decoherence induced by an environment 
containing a few discrete modes has received attention 
in recent years in the study of quantum information pro- 
cessing using wave packets. There are many situations 
where the environment has a discrete spectrum. The 
decoherence of molecule vibrational wave packet due to 
the vibrational-rotational coupling is a typical example 
[H [l^]. Other physical situations where environments 
have only a few modes being relevant include electromag- 
netic field in cavity QED, molecular qubit decoherence in 
glassy environments [2l| . and electron transfer and exci- 
ton decoherence dynamics in proteins [2^ , to name a few. 
In fact, many environments associating with manipula- 
tions and measurements in quantum information process- 
ing involve a finite and discrete modes. The decoherence 
process under such circumstances must be studied with 
the back-reaction of the environment to the system be- 
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ing fully taken into account. The result is expected to be 
different from that happens with a continuous spectrum 
environment. 

For making comparison between the mechanisms of de- 
coherence in a few discrete modes and a continuous spec- 
trum environment, we take the Caldeira-Leggett model 
Q with a discrete spectrum environment as a representa- 
tive one. In this model, the principal system is a particle 
in a harmonic trap and the environment is a set of nonin- 
teracting harmonic oscillators each of which linearly cou- 
pled to the principal system in their coordinates. Many 
people have studied this model, some have done research 
on models slightly modified from this concerning topics 
not only of decoherence but many other issues of o pen 
quantum systemsjlji, S EE IHIil, El, Q M 

H [O, E HI [mISI [lElU An exact master equa- 
tion to the motion of the central oscillator with a general 
spectral density at an arbitrary ternperature was given 
in [Toll and derived alternatively in p8| . The exact mas- 
ter equation lies on the dissipation dynamical equation 
which fully takes into account the back-action effects of 
the environment on the system. Here we present an an- 
alytical solution to the dissipation dynamical equation 
to directly describe the non-Markovian processes of the 
central particle. And the wave packet dynamics can be 
expressed in terms of the solution to the dissipation dy- 
namical equation without invoking the master equation. 
As a result, a qualitative change of the dynamics of the 
central particle is directly read from the exact solution 
to the dissipation dynamical equation we obtained. The 
resulting effect can be used to analyze decoherence and 
instability of the system under the influence of the envi- 
ronment. 

The paper is organized as follows. In the next section, 
we briefly describe the Feynman- Vernon influence func- 
tional approach to a quantum harmonic oscillator inter- 
acting linearly with a general environment of N modes. 
We then in the following section solve analytically the 
reduced dynamics of the system. We obtain two intrinsic 
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energy scales relating to the time scales of the system 
and the environment. We will explicitly show how the 
different energy scales between the system and the envi- 
ronment alters qualitatively the dynamics of the system. 
We also take the continuum limit of the environmental 
frequencies to the Ohmic spectrum |23l| and recover the 
solution previously obtained in 0, In section IV we 
study the evolution of a single wave packet. We introduce 
a decoherence measure for the single wave packet evolu- 
tion. Non-Markovian decoherence of the wave packet in 
such an environment with a few discrete modes is demon- 
strated by a modified Drude's spectral density ,29]. Phys- 
ical implications of our results will be given and discussed 
in the conclusion section. The detailed mathematics is 
left in the Appendix. 

II. TIME EVOLUTION OF OPEN QUANTUM 
SYSTEMS 

Following many works on Caldeira-Leggett model [3, 
0] , we briefly review the main results obtained previously 
that will be used later in this work. The Hamiltonian of 
a particle in a harmonic trap linearly coupled with an 
environment is given by 

N 

+ Y.C,Xq,. (1) 

The first two terms in the bracket is the Hamiltonian 
of the system and the second term as a summation is 
that of the environment consisting of a finite N discrete 
modes. The third term is the interaction Hamiltonian 
between the system and the environment. The notations 
follow the convention. The requirement of a finite N 
discrete modes in the environment makes our consider- 
ation different from original Caldeira-Leggett model. It 
should be particularly noticed that in Caldeira-Leggett 

model an additional counter-term, „'^' a X'^, is usu- 

ally added to Eq. ^ in order to study the generic be- 

I 



where 



N p2 

Kl{T-T') = -y^^P^s\nUOk{T-T') (6) 



havior of dissipation dynamics for a harmonic oscilla- 
tor. Such a counter-term exactly cancels the coupling- 
induced frequency-renormalization effect 0. Since our 
interest lies on the environment with a finite discrete 
modes where no high frequency cut-off is introduced so 
that no counter-term is needed for Eq. |[T]), according to 
the standard renormalization theory 30]. As it has also 
been pointed out [6i, ,31] the coupling-induced frequency- 
renormalization effect in many cases is a physical observ- 
able effect that has to be taken into account. In fact, the 
study of quantum decoherence should fully explore the 
environment induced dissipation, fluctuation as well as 
renormalization effects on the system. Nonexistence of 
the counter-term for the present case of a discrete spec- 
trum environment enables us to explore the decoherence 
dynamics of wave packets purely induced by the interac- 
tion with the environment without ambiguity. 

By using the common assumption that the system and 
the environment is initially uncorrelated while the en- 
vironment starts with an equilibrium state, the reduced 
density matrix of the system evaluated in the position 
basis at time t is given by [1, 0] 

< x\PA{t)W >= PA{x,x',t) 

= j dxodx'aJ{x, x' , t\xo, x'^, 0)pAixo, xf,, 0), (2) 

where J{x, x' , t\xQ, x'q, 0) is the propagator which includes 
the back-action effect of the environment on the system 
and is given by 

J{x,x',t\xo,x'o,0) = yP[x]P[x']e^(^^W-^-^M).F[x,x'], 

(3) 

in which Sa is the action of the system 

and J-[x, x'] is the influence functional describing the in- 
fluence of the environment on the system. The explicit 
form of the influence functional has been well solved i2J| : 



(5) 

I 

and 

^ii(T-r')-E:r^coth^cos^fc(r-r'), (7) 
^ 2mwfc 2 

are defined as the dissipation and fluctuation Kernels re- 
spectively, and P is reciprocal of the product of the initial 
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equilibrium temperature and Boltzmann constant. 

Eqs. (I5][5]) shows that the cfFective action has a 
quadratic form. Hence the path integral can be exactly 
carried out with the stationary path method . By in- 
troducing the new variables R(t) = ^(^)+^ ('^) ^^(y) = 

I 



where Xijit) ihj — li2) is defined by 

X^At) = f dT r dT'v,{T)KR{T - t')Vj{t') (9) 

Jo Jo 

and Ui and Vi satisfy the following dissipation dynamical 
equations of motion, 

u,(r)+r!V(r) + -^ / dr'if/(r~T')u»(r') = 0, (10) 

i),(T) + n^v^ir) - ^ ^ dr'Kjir - r')t'.(r') = 0,(11) 

subject to the boundary conditions wi(0) = ""1(0) = 

V2{t) = U2{t) = l,«i(t) = lii(t) = V2{0) = U2{0) = 0. N„ 

is the contribution of fluctuations around the stationary 
pathes which is independent of the position variables for 
such a quadratic form of the Hamiltonian. 

III. ANALYTICAL SOLUTION TO THE 
DISSIPATION DYNAMICS 

Since Eq. (jlip is the backward version of the forward- 
ing equation Eq. (flO)) . the propagator of the system, 
Eq. ([8]), is completely determined by solving the dis- 
sipation dynamical equation (|10p . The first two terms 
in Eq. pU)) correspond to the classical equation of mo- 
tion for the coordinate of the central oscillator and the 
third term with the dissipation kernel Kj accounts for 
the non-Markovian (memory) effect of the back-action of 
the environment on the system. While, the fluctuation 
contribution from the environment is given by Eq. Q. 

On the other hand, the propagator Eq. ([8]) is built upon 
the solutions of the dissipation dynamical equations (fTUl 
[TT|) . The solution of these equations gives rise to the full 
dynamics of the system. Using the fact that Eq. (fTT|) 
is the backward version of Eq. (jlOp , vi (r) and V2 (t) are 
related to ui(t) and U2{t) by wi(t) =^ U2{t — t), V2{t) — 
uiit — r). Furthermore, ui{t) and U2{t) can be expressed 
as mi(t) = Z{t) — ^^Z{t) and U2{t) = where Z{t) 
obeys the same equation as Eq. pH)) . 

Z{t) + n'ZiT) + A ^ dr'KiiT - t')Z{t') = 0, (12) 



x(t) — y (r), we rewrite the propagator J{x, x' , t\xQ, Xq, 0) 
as JCR, r, t|i?o, ^'o, 0). The resulting propagator becomes 



(8) 

I 

with the boundary conditions ^(0) = and 2'(0) — 
1. The solution of Z{t) determines the behaviors of 
wi(t), U2(j) and V\{t), V2{t) and all the relevant dynam- 
ics of the reduced system. Z{t) was solved in cases of 
continuous spectra of Ohmic and super-Ohmic types at 
the Markovian limit 33] . Here we shall give the analyti- 
cal and exact solution of Z{t) for discrete spectra. 

We found that the feature of Z{t) depends on a rela- 
tionship between two intrinsic energy scales: MQ,^ and 

EH ^ The former is the intrinsic energy scale of 

the system in terms of the mass M of the central parti- 
cle and the corresponding frequency $1 in the harmonic 
trap. The latter is associated with an interaction energy 
scale between the system and the environment character- 
ized by the coupling constants {Ck] k — 1, • • • , N} and 
the intrinsic energy scales of the environmental oscilla- 
tors {mLol]k = 1, ■ • ■ ,N}. We call Mfl'^ the bound- 
ing strength of the system because it is just the second 
derivative of the harmonic potential at its valley. This 
strength tells how strongly the particle is bounded in this 

harmonic trap. We name X^feLi simply the spectral 
strength. These two energy scales indeed act as two time 
scales in the total system, one corresponds to the time 
scale of the system and the other corresponds to the time 
scale of the environment accompanied with the couplings 
between the system and the environment. Since the en- 
vironment contains many different frequencies, it is not 
easy to figure out a definite time scale. However, the 

spectral strength ^^IJ? contains the contribution of 

every individual frequency in the environment plus the 
corresponding coupling to the system. It naturally pro- 
vides an alternative expression to the time scale of the 
environment. 

We will see that different relationship between these 
two strengthes (or equivalently between two time scales 
of the system and environment) alters the overall form 
of the solution to the dynamics of the system. The ex- 
plicit forms of Z{t) in cases that the bounding strength 
of the system is stronger than, weaker than or equal to 
the spectral strength are quite different. We just present 
here the results and leave the verification of the solutions 
Z{t) in Appendix. In general, Z{t) has the following 



J{R,r,t\RQ,ro,0) = f^^f,{TiMlroB.oU2it)-roB.U2iO)+B.oruiit)~rRui{0)]+ilxil(t)rf, + ^^^^ 
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form 



N 



fc=0 11 



n,=i('^. 



J'fc) sin (i/fcr) 



(13) 



where z/„ satisfy 



N 



2—1 k—1 i^k.i—1 

(14) 

forn = 0, • • ■ , iV. When Mn'^ > J2k=i 



^ we have 



(15) 

For a discrete spectrum, this constitutes an oscillation 
function in time. The continuous limit will lead this 
summation over sine functions to a Fourier integral that 

may not always be oscillatory. When MV,^ — X^fcLi 
— and the solution becomes 

N 9 



nuJf 



-1 ''k 



k=l 



N 



rijliK'-J^fc) sin(i/fcT) 



(16) 



with all the other f„'s, iV > n > 1, still being located 



by (Uni). In the case that Mn^ < J2k=^ ^ 
where /iq is a real number and Z(s) becomes 



1^0 = «Mo 



nf=iK+A*o) sinh(AioT) 



Z<ir) 



n^r(-: 



N 

E 



2 

i 



sin (vkT) 



1 K+/^o)njlij^fe('^: 



(17) 



As we have seen when the bounding strength of the 
system is stronger than the spectral strength, the oscil- 
latory nature of the system is maintained although its 
oscillation details are completely altered. Coupling to 
the N oscillating modes of the environment generates N 
new oscillating modes to the system's dynamics. Com- 
bining with the original oscillating mode f2, the system 
has now totally -I- 1 oscillating modes, as shown in 
Eq. p3)) . Note that the existence of -I- 1 new oscil- 
lating modes was also obtained by Haake and Reinhold 
through a normal mode transformation to the model 
Hamiltonian but the physical picture is different. In Q, 
the A^ + 1 new oscillating modes is just a hybridization 
to the original frequencies {ft, ujk, k — 1, ■ ' ' i ^} of the 
A" -|- 1 harmonic oscillators for the system plus the en- 
vironment, due to the coupling between them. Here the 
solution of the system, after the environment's degrees 
of freedom are integrated out, contains A^ -I- 1 oscillat- 
ing modes. It indicates that after the environment is 
traced over, the system is no longer a particle trapped 



in a harmonic potential, as manifested by the solution of 
Eq. (Uni). The N+1 oscillating modes n = 0, • • • , TV} 
are determined by Eq. p^ . 

Furthermore, when Mfi^ y^t^i we find from 

Eq. (fTil) that the frequency vq 0. In particular, the 
condition Mfl^ = X^feLi ^^JJ" leads exactly to vq — 0. As 
we see, the first term in Eq. (|16p comes from i^o = 0, while 
the rest A^ mode oscillations in the solution of Z{t) are 
induced by the coupling of the system to the A^ modes 
of the environment. This solution indicates that at the 
critical energy condition, Mfl^ = Y^^' , the central 

particle is driven out of the harmonic potential by its 
interaction with the environment. After a relative long 
time, the solution of Z{t) is dominated by the linear term 
in time, which is responsible for an irreversible process 
as a sign for the possible rise of forever loss of quantum 
coherence. 

Now, we further consider the case of Mft^ < 
J2k=i interesting to see from Eq. (fT4|) that vq 

becomes an imaginary number. Although the solution 
Eq. (fT7|) still contains A^ oscillating modes induced from 
the coupling of the system to the A^ modes in the environ- 
ment, the imaginary root — ifiQ provides the solution 
of Z{t) with a component that exponentially grows up 
in time. The stronger the spectral strength is, the larger 
/io will be (see the appendix A.l). In other words, when 
the bounding strength of the system, MV.^, is below a 
critical value of the environment's spectral strength, the 
system will be pulled out from the harmonic bounding 
potential very quickly and its dynamical process becomes 
irreversible in an exponential growth rate. 

In fact, the existence of two energy scales in the Hamil- 
tonian Eq. ([T]) and the corresponding non-Markovian dy- 
namics resulting from the competition between the two 
energy scales have not be paid attention in the literature. 
This is mainly because a counter-term, Yli^i „'^' ^ X^, is 

usually added to Eq. ([1]) in order to study the generic 
behavior of the dissipation dynamics for a damping har- 
monic oscillator Q , or a " positivity condition" , 17^ — 

X]i=i ^ ^ (for the case M = m), is imposed on 
the coupling constants and the unperturbed frequencies 
to ensure the Hamiltonian having a finite lower bound 
Q- Obviously the counter-term or the positivity con- 
dition excludes the dynamics corresponding to the case 



l^k= 



However, in the present work we 



deal with the environment having only a finite number of 
discrete modes in which no high frequency cut-off needs 
to be introduced. Therefore, no counter-term exists ac- 
cording to the standard renormalization theory [3d| . Also, 
as Caldeira and Leggett had extensively discussed 
the coupling-induced frequency-renormalization effect in 
many cases is indeed a real physical effect that has to be 
taken into account. There are many physical situations, 
such as phase transitions of early universe in extreme en- 
vironments in cosmology |37] , fission of a heavy nucleus 
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with coupling of the collective degree of freedom to the 
single-particle modes in nuclear physics 0, atomic tun- 
neling with phonon coupling in glasses [33, and the in- 
stability of Bose-Einstein condensation influenced by the 
trapping field in cold atoms ji9j, etc. where the renor- 
malization effect can be very large to reach the regime 
Mfi^ < Vf' , such that the system can render the 
original potential minimum unstable. The corresponding 
dynamics will undergo a catastrophic change that can be 
observed experimentally. Therefore the positivity condi- 
tion should be also not applicable to these physical situ- 
ations, in contrast to the study of environment-induced 
dissipative or damping harmonic oscillation in the liter- 
ature 0,[2i|. 

In order to make a comparison between the environ- 
ments with discrete and continuous spectra, we take 
the continuum limit of the spectrum of the environ- 
ment. We first define a polynomial whose roots are 
those frequencies i^k, k — 0, ••• ,iV. This is the poly- 
nomial of Eq. m]), P{iy) = Mp(v)W{v) where p{v) = 

We observe that the general form of Z{t) can also be 
rewritten as Z{t) = X]fc=o p^(uk) ^^^(^fc''")' where p' de- 
notes the derivative of p with respect to v. By the 
fact that the roots of P{v) are also that of p{v) and 
the roots are paired due the evenness of the poly- 
nomials, we further rewrite Z{t) replacing sinivj^r) 
by — — as Z(t) = — i '5^?.;'^^"'^ ttttA^ enumerat- 



'^f^ as Z{t) = 

ing —Vq — = l'N+2, - ■ ■ ,^t^N = J^N+N+l- 

By the residue's theorem, Z{t) becomes ^ dz '^^^^ 
where the integration contour encloses all the roots of 
p{i'). Then taking the continuum limit of p{i') that 



N Cf 



D{uj)C^{uj) 



where D(uj) 



J = l m uj'j—u^ ' Jo 

is the density of states in the environment and setting 
= 2^70^9 (A - uj) as the spectral density of 
an Ohmic environment (here Q is the Heavyside func- 
tion), we reproduce 



Z(t) 



— p-70 



^ sin(j7T) 



(18) 



as a solution to the corresponding version of dissipation 
dynamical equation Eq. Z{T) + 2joZ{T)+nlZ{T) = 
0. Here 70 is a constant that is usually very small in 
comparison to the high frequency cutoff A, and f2 = 



7o 



fir (for 7o <C A) with the system's reno- 



mahzed frequency fl^ = fl^ - Ajo^/tt 0, US, El. As 
we see, whenever the environment's spectrum becomes 
continuous, the form of Z{t) as a summation over sine 
functions in Eq. ([13]) becomes a Fourier integral and it 
may not always be oscillatory. In other words, the inter- 
action of the system with a discrete spectral environment 
leads to very different results from the system interacting 
with a continuous spectrum environment. 

Having obtained the analytical solution to the dissipa- 
tion dynamical equation for an environment with discrete 
spectrum and reproduced the previous known results for 
the Ohmic bath in the continuous spectrum limit, we 
can now utilize this analytical solution to study the wave 
packet dynamics in the next section. 



IV. NON-MARKOVIAN WAVE PACKET 
DYNAMICS 



The geometry of a wave packet is well described by the 
covariation matrix [3J| : 



AX^ A{XP} 
A{PX} AP2 



where AX'^ and AP^ are the widthes of the wave packet 
in position and momentum spaces respectively, and are 
defined via AX'^ = (X^) - (xf, X^ = (X), AP^ = 

(p2^ - {Pf and Pc = (P) with the bracket (.) denot- 
ing the expectation value of an operator in the state pA- 
A{A"P} is a quantity that measures the internal corre- 
lation between the position and momentum observables, 
A{A:P} = ^{XP + PX) - {X){P) = A{PA:}. The 
density matrix of the initial wave packet of the system 
generally takes the form: 



fU n\ AT -J^[iRo-iA{XP}orof + -^AX^AP^r^„+2XoRo + ^{Po~A{XP}„Xo)ro\ 

PA[Ro,ro,0) = Noe ^^^0 '- 



(19) 



r 



where iVo is the normalization constant, and the sub- 
script " 0" here denotes the initial values of these quanti- 
ties. 



The density matrix at time t is determined by the time 
evolution equation, 

PA{R,r,t) ^ / dRodroJ{R,r,t\Ro,ro,0)pA{Ro,ro,0) 



(20) 

where the propagator J{R,r,t\Ro,rQ,0) (see Eq. ([5])) is 
easy to obtain once the dynamical equations (jlOlllip are 
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solved exactly. Explicitly, the solution of Eq. (1201) 



^[{B~j-A{XP}it)r)\-^{AX^(t)AP^(t))r^+2X,{t)R+S^{p4t)-A{XP}(t)X4t))r] 



(21) 



r 



which keeps the same form as the initial wave packet ([191) > 
and iV is a normalization factor. The time evolution of 
the wave packet is thus completely determined by the 
time evolution of the covariation matrix elements, which 
can be found from Eqs. ^ and (HOj)- Together with the 
solutions of Z{t) discussed in the last section, the full 
information of the wave packet dynamics, in particular, 
the decoherence process can be analyzed now. 

The decoherence behavior is mainly characterized by 
the decay of the off-diagonal element of the reduced den- 
sity matrix. In the literature, to make the quantum co- 
herence of wave packets manifestation, one usually starts 
with a superposition of two well-separated Gaussian wave 
packets where the decoherence can be measured dir ectly 
from the interference of two wave packets 0, [H, Il4| . 
However, a single wave packet itself represents a macro- 
scopic quantum state which has its own physical interests 
in many physical systems described by wave packets. A 
decoherence measure to a single wave packet is desirable. 

Recall that pA{R,r,t) = {x\pA{t)\x') with R = {x + 
x')/2 and r = x — x' , the off-diagonal matrix element 
is given by the r-dependent part in (|2ip which describes 
the quantum coherence dynamics of the wave packet. In 
fact, Eq. (pil) tells that the dynamics of the off-diagonal 
matrix element is fully determined by the covariation 
matrix AX'^{t), AP^{t), A{XP}{t). In particular, the 
correlator A{XP}{t) describes the phase dynamics of 
the density matrix, while the Heisenberg uncertainty, 
AX'^{t)AP'^{t), measures the amplitude dynamics of its 
off-diagonal behavior. Note that because of its inde- 
pendence on the widths AX'^{t) and AP^{t) as well as 
the correlator A{XP}{t), the center motion of the wave 
packet (given by Xc and Pc) does not change the geom- 
etry of the wave packet. Without loss of the generality, 
we may let the initial position and the momentum of the 
wave packet be zero, Xq ~ Pq = 0, then the reduced 
density matrix simply becomes: 



PAiR,r,t) = Ncxp 



R' 



A{XP}{t) 



2AX^{t) h AX^{t) 
/ AX\t)AP^t) - AjXPyjt )^ 

V WKxHt) 



Rr 



(22) 



Now it becomes clear that the decay of the off-diagonal 
matrix element is determined by the quantity D^it) = 
AX (t)AP^(^-A{xp} (t) ^ ^Yiere the denominator AX^{t) 
measures the delocalization of the waver packet in- 
fluenced by the environment, while the numerator 
AX^{t)AP^{t) ~ A{XP}^{t) is actually the Schrodinger 



uncertainty [35l|. The minimum Schrodinger uncertainty 
has been used as a criterion to examine if a wave packet 
is a squeezed coherent state ^3]. Thus Dc{t) is a natu- 
ral quantity to characterize the delocalization of a wave 
packet accompanied with wave packet decoherence. The 
degree of quantum decoherence in a wave packet can then 
be extracted from the off-diagonal matrix element: 



1 AX^(t)AP-'(t)-A{XP}^(t) „2 



(23) 



The larger Dc{t) is, the less there is quantum coherence. 

To demonstrate the non-Markovian decoherence dy- 
namics of the wave packet, we plot the time evolutions 
of the wave packet width and the corresponding deco- 
herence measure at different physical conditions. Note 
that for an environment with a few modes, the parame- 
ter (3 bares no sense of thermodynamics and it is just a 
parameter in the initial state of the environment. We 
set here AI = m = n = h=l, (3~^ = 1.15 and 



Ck = 



r2 



with r = 500. The param- 



os 



eter 7 is dimensionless and is set differently according 
to the relationship between Mfi^ and Vf^ , We 

take an initial wave packet as AXq — 1/5, APq = 5 
and A{XP}o ^ Xq ^ Po ^ 0. We first fix the fre- 
quency distribution in the environment and vary only 
7 to obtain the conditions Mfi^ > ^ 



Z^fe = l 



-2% and Mfl^ < Vf 1 We then fix 7 and 

change the frequency distributions to satisfy the above 
three relationships between the bounding and the spec- 
tral strengthes. The results are shown Figs. 1-2. From 
the plots, one can see that when the bounding strength 
is larger than the spectral strength, the wave packet 
keeps oscillating in its width as well as the decoherence 
measure. There is no unidirectional growth of the de- 
coherence measure. When the two strengthes balance, 
the wave packet starts to spread. But the decherence 
measure still oscillates at some finite values. The mono- 
tonic loss of quantum coherence is seen when the spectral 
strength wins over the bounding strength. This happens 
also together with an even faster spreading of the wave 
packet. These phenomena are seen no matter 7 is var- 
ied with the frequency distribution fixed or the frequency 
distribution is changed with 7 fixed (see Fig. 1, and the 
solid-curves in Fig. 2). Consequently, the emergence of 
decoherence is subjected to the definite greatness of the 
spectral strength in comparison to the bounding strength 
of the system here. If the spectral strength is just equal 
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to the bounding strength, delocahzation dynamics occurs without necessarily leading to decoherence. 
I 
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FIG. 1; (a), (b) and (c) show the time evolution of the width AX^ of the wave packet at the condition AIQ^ >,= and 

< '^Z^y respectively, (d), (e) and (f) the time evolution of the decoherence measure Dc at the corresponding conditions. 

The frequencies are chosen as {a;i, W2, ws, u;4, a;5}={0. 48, 0.86, 1.72, 1.84, 1.89} . The coupling constant 7 is adjusted as 0.01, 
0.39 and 0.58 to satisfy the bounding strength being larger than, equal to and smaller than the spectral strength. The wave 
packet delocalizes whenever its bounding strength is no larger than the spectral strength. Dc grows up monotonically when 



k — 1 rnuj'f 



The balance condition AID, 



EN 
fe=l 



results in slow diffusion of the wave packet, slow compared to 



that for Mfi^ < Y.k^i 



ct 



without rendering unidirectional loss of quantum coherence. 



The non-Markovian decoherence dynamics of the wave 
packet (the memory effect induced from the back-action 
of the environment) can be seen by contrasting the exact 
solution with Markovian approximation which is phys- 
ically valid when the time scale of the environment is 
much smaller than that of the system and has been 
widely used in the literature [2^, H^]- Mathematically, 
the Markovian approximation takes the history indepen- 
dence of the dynamical variable in Eq. ()12|) . namely 

/J" dT'Kl{T-T')Z{T') ~ Z{t) /J' dT'Kj{T-T') lU Solviug 

the dissipation dynamical equation [3G]. When the time 
scales of the environments and the system are compara- 
ble (equivalently when the values of the two strengthes 
approach to each other) , we expect the dynamics of the 
wave packet under Markovian approximation deviates 
obviously from the exact solution. The three frequency 



distributions are chosen such that the time scales of the 
environments fall upon the regimes of interests. 

As we can see from Fig. 2, when the environment con- 
tains only frequencies much larger than the central one, 
the system is located in the regime AIQ^ ^ ^k=i ^^J^"- 
We find that the Markovian approximation well agrees 
with the exact solution. However, when the frequen- 
cies of the environmental oscillators all lie near the cen- 
tral frequency (strongest memory regime) where the two 
strengthes are very close to each other, an apparent dif- 
ference between the Markovian and non-Markovian re- 
sults shows up. The exact time evolution of the wave 
packet width exhibits several kinks and bumps that are 
smeared when Markovian processes are asserted. The 
extension of wave packet delocahzation is over estimated 
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by Markovian approximation. The exact time evolution 
of the decoherence measure in this time scale shows os- 
cillatory features that are not seen from the Markovian 
approximation. The magnitude of the decoherence mea- 
sure is also over evaluated in the Markovian approxima- 
tion. In other words, ignoring the history dependence 
smears out the subsequent features of the time evolu- 
tions of AX^ and in this case. When it only consists 



of lower frequencies, namely the environment has a longer 
time scale than that of the system, the system is easy to 

fall into the regime Mfi^ < Vf!' , where we find the 
Markovian approximation over evaluates the magnitudes 
of AX'^ and Dc but not as severe as that in the balance 
case. 




The wave packet dynamics we obtained here is not 
limited to environments with discrete spectra. We may 
consider the wave packet evolution in an Ohmic bath (a 
continuous spectrum environment) as an example. In 
this case, the system undergoes a underdamping process 
according to We find that when t ^ oo the width 



of the wave packet AX'^{t) is given by 

AX^t ^ oo) = - / dwcoth(/3— ) 
Jq 2 

M ] (24) 

Vm (r!2-u;2)2+ 472^2 ^^'^^ 

Note that in the linear response theory, a harmonic 
oscillator of frequency fir reacts to an external force 
with the imaginary part of its response being <t{uj) = 
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{jT (n^^-Jy+A-y-^u^ ) ■ ^ rewritten as 

AX2(oo) ^ I /~dwcoth(/3%i )cr(tj) which is a well- 
known result of the fluctuation-dissipation theorem 0, 
|40| . In the extremely underdamped limit, the spread 
of the wave packet ultimately goes to an equilibrium 
value: AX'^{oo) = jmr coth(/3^), the same result ob- 
tained in the equilibrium quantum statistical mechanics 
SH- The numerical evaluations of AX'^(t) at large 
times under such a circumstance also accurately match 
the equilibrium quantum statistical result. 



IV CONCLUSIONS AND DISCUSSIONS 



quent wave packet dynamics. The over all trends of wave 
packet dynamics are not altered by Markovian approxi- 
mation when the bounding strength is much larger than 
the spectral strength. But Markovian approximation un- 
der the condition of a larger spectral strength changes 
the wave packet dynamics very drastically from the exact 
solution. The critical condition that the two strengthes 
balance each other manifests the Non- Markovian wave 
packet dynamics most significantly. In conclusion, we 
have taken a few factors to explore the mechanisms of de- 
coherence of a wave packet in a few-mode environment. 
One may attempt to see these issues in relevant systems 
mentioned above. 



Here we summarize what we have done in this paper. 
We study the dynamics of a particle in a harmonic trap 
linearly coupled to a set of noninteracting harmonic os- 
cillators. We obtained an exact analytical solution of 
the system with the environment containing discrete N 
modes. By taking the continuous spectrum limit, we ex- 
amined the dynamics in an Ohmic spectral density and 
reproduced the solution once obtained by others. In order 
to visualize the dynamics more vividly, we let the system 
start with a wave packet. The decoherence processes are 
described in terms of the geometry of the wave packet. 
We then closely looked at the wave packet's localization 
and delocalization behavior against its decoherence mea- 
sure defined by the ratio of Schrodinger uncertainty over 
the delocalization extension of the wave packet. We ob- 
serve two intrinsic energy scales, the bounding strength 

MV,'^ and the spectral strength X^fcLi acting as the 

time scales of the system and the environment, respec- 
tively. These two energy scales are veiled in the previ- 
ous investigations when a counter-term was added to the 
Hamiltonian in the continuous spectral density where a 
high-frequency cut-off has to be introduced or a positiv- 
ity condition is imposed to the system. However, the 
relationship between these two intrinsic energy scales de- 
termines the feature of the delocalization as well as deco- 
herence dynamics of the wave packet in the cases where 
the environment has only a few modes so that no counter- 
term exists or no positivity condition can be added when 
decoherence dynamics is concerned. 

We also numerically demonstrated the time evolution 
of the wave packet based on the exact solutions we ob- 
tained with an environment of a few modes satisfying 
a modified Drude's spectral density. Our results show 
that a stronger spectral strength drives the wave packet 
out of the harmonic bounding potential as a consequence 
of the environment-induced instability, which results in 
strong entanglement between the system and the envi- 
ronment and leads to severe decoherence, while a mod- 
erate or weaker spectral strength oscillates its decoher- 
ence measure over time. We also compared the exact 
results with Markovian approximation. If Markovian ap- 
proximation is done ahead, the relationship between the 
two strengthes still plays a determinant role in the conse- 



Acknowledgments 

We would like to thank Drs. M. T. Lee and J. H. An 
for fruitful discussions. We would also like to specially 
thank Prof. A. J. Leggett for his help in clarifying the 
counter-term problem in Caldeira-Leggett model. This 
work is supported by the National Science Council of 
Republic of China under Contract No. NSC-96-2112-M- 
006-011-MY3 and No. NSC-95-2112-M-006-001. 



APPENDIX: THE SOLUTIONS TO THE 
DISSIPATION DYNAMICAL EQUATIONS 

A.l Root Property 

Before we proceed to justify our solution to the dis- 
sipation dynamical equation, we first look at (jl4p . The 
roots of this polynomial determine the properties of Z(s). 
We denote the general argument of this polynomial by v 
instead of i>n- For convenience of reference, we write this 
polynomial again, 

^oM = 

N N N 

i—1 k—1 i^k,i—l 

(25) 

It is well known that if the polynomial is evaluated at 
some point on the real axis to be positive and is negative 
at some other point, then there is a point in between 
at which this polynomial is evaluated to be zero as its 
root, denoted by Vi. We arrange lui through wat in a way 
that < Li;2 < • • • < Wat and let v = ujn, then Eq. (^51 
becomes 

PoK)=(-ir— n(^n-^') n (^f-^^)-(26) 

i—l i—n-\-l 

Clearly the sign of the polynomial at cun for n = 
1, 2, • • • ,N is determined by (— l)", so the roots are lo- 
cated in the vicinity between adjacent frequencies Un and 



uJn+1 and we have located N — 1 distinct roots that are all 
positive. Since the polynomial is even in its argument, 
the roots are paired. This gives another set of — 1 

I 



Po{{aN+m) = 



We can make Po{aNil, = un) and Po((aw + f)^) differ 
by a sign by choosing an / which is large enough that, 

and the sign of Po((aw + is (-1)^+^ So far by 

including the negative partner of this root, we have found 
2N distinct roots. We know this polynomial has + 2 
roots so there is still a pair of roots not being located yet. 
Since the sign of the polynomial at lo\ is minus, if -Po(O) 
is positive, then there is a root lying between and uj\. 
We see -Po(O) is just the difference between the bounding 
strength of the system and the spectral strength times a 
positive number. 

If MVl^ is larger than then Po(0) is positive 

and we have a root that is larger than and smaller than 
Since -Po(O) is proportional to the product of all the 
roots it has, Po(0) = implies there is a pair of roots 
doubly located at the origin. If Po(0) is negative, then 
we have to look for roots somewhere out of the real axis. 
In this case we define Py{y) = Poiw) where y is in the 
real domain, then Py{0) — Po(0) < 0. If we can find a y 
such that Py{y) > 0, then there is a root between and 
y for Py which implies a root for Pq. Let y = fil where 
/ is a positive number and we have 

N 

Py{m^l[[{na,f + ifnr)]if + i) 

gi;r (Oo,)^ + (/n)0 '^^) 
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roots that are all negative. We can also see there is a 
root larger than wjv by letting uji — aiQ for i = 1, ■ ■ ■ , N 
and v = (ajv + f)^, that is. 



I 

Again, / can be chosen to be large enough such that 



This gives us a real root to Py and also an imaginary root 
to Pq. This root of course has a negative partner and we 
find finally all the roots. 

We summarize the location of half of the roots as 



< i^a < uji < < UJ2 < ■ ■ ■ < i^N-i < < i^N < oo, 

(30) 

for Mfl^ > X^^i and vq — z/zq where /xq is a real 

number for Mfl^ < X^A^i The other half is ncga- 

tively paired with this. 



A. 2 Verification of the Solution 



We have proved in the last subsection that j/q through 
are + 1 distinct numbers. This ensures that the 
amplitudes of each component in the solution of Z{s) 
do not diverge. Now we shall check whether fits 
p2| with the boundary conditions that Z{0) — and 
Z(0) = 1. Substituting ^ to we have 



(-1)^+1 lliiaj, + If - a]){{a^ + jf {m^^ 



•^N C 



1 



(27) 
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(Z{s) + n'Zi.s)) + - / ds'Kj{s-s')Zis') = 

I n 



fc=0 
- 



N 



- sin(cjjs) 



(3i; 



where AvS are the aniphtudes — ] i'' We 
now introduce an identity, 



E 



r 0, 1 < j < ft - 2 

= i 1, j = n-l (32) 



Comparins; the factor Vf^ „ l^'w — ^ in the last 

hne of PT|) with the identity ([5^ by replacing 6^ by vf^-^ 
and n by + 1, we can see in the numerator that the 
highest power oi i^l is N - I and A^-l = (iV + l)-2== 



n — 2. So this summation goes to zero. And Z{s) with 
the proposed form (fT3)) fits (fT2|) . The boundary condition 
Z(0) = is easy to check for sine functions are zero when 
their arguments are zero. To check the boundary value 
satisfied by its first derivative, we take, 



N 



Z{0) = Aki^k = V 



(33) 



Making use of ([32|) for the highest power in the numerator 
now equaling to A'^ n — 1, Z{Q) becomes unity. 
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